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We study an anyon model in a toric honeycomb lattice. The ground states and the low-lying 
excitations coincide with those of Kitaev toric code model and then the excitations obey mutual 
semionic statistics. This model is helpful to understand the toric code of anyons in a more symmetric 
way. On the other hand, there is a direct relation between this toric honeycomb model and a 
boundary coupled Ising chain array in a square lattice via Jordan- Wigner transformation. We 
discuss the equivalence between these two models in the low-lying sector and realize these anyon 
excitations in a conventional fermion system. 
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I. INTRODUCTIONS 

Kitaev toric code model opens a new direction to study 
the storage of quantum information and design of quan- 
tum computation in a topological protected way [l[ . The 
objects carrying the quantum information are so-called 
anyons, the quasiparticles obeying exotic statistics in 
two-dimensional condensed matter systems [I, Q . 

The first realistic anyonic quasiparticle is Laughlin 
quasiparticle in v = 1/S fractional quantum Hall sys- 
tem and possibly obey 9 = fractional statistics 
0, Q- Nowadays, the researches for anyons mainly focus 
on the following aspects: non-abelian fractional quantum 
Hall states |9l4lll|: topological insulator carrying Majo- 
rana fermionic edge mode [l2l |: and Kitaev lattice spin 
models, i.e., the toric code model Levin- Wen model 
[l3j and honeycomb- lattice spin model [l^ . 

The abelian anyons in Kitaev models have currently 
attracted many research interests. Experimentally ex- 
citing, manipulating and detecting abelian anyons have 
been suggested or tried for the toric code model [l5j 
and for the insulating phase of Kitaev honeycomb-lattice 
model [l6j . The explicit presentation of nonabelian 
anyons has also been studied [T7| . 

The realization of Kitaev toric code and honeycomb 
lattice models is not easy because of the unconventional 
interaction between spins. In this paper, we propose a 
toric honeycomb lattice model which is equivalent to Ki- 
taev toric code but in a more symmetric way. It was 
known that the insulator phase of Kitaev honeycomb lat- 
tice model is equivalent to Kitaev toric code model [l^ 
according to a perturbation analysis. 

Instead of such a complicated way, in this paper, we di- 
rectly use the group element of Z2 gauge group in Kitaev 
honeycomb model to construct our model Hamiltonian. 
There are two independent operators to represent the Z2 
group. They play roles of the stabilizer operators of the 
toric code. Thus, we find that our model is equivalent 
to Kitaev toric code model but the Hamiltonian is more 
symmetric. All the eigenstates can be known because 
the model is exactly solvable. The low-lying excitations 



contain a closed subset in which the excitations obey the 
mutual semionic statistics. 

Furthermore, we show that this toric honeycomb lat- 
tice model can be mapped to a decoupled Ising chain 
array in a square lattice under a special external field. 
The latter has been studied before by one of us with 
Li [1^ and the low-lying excitations of that model also 
obey mutual semionic statistics. However, the decoupled 
Ising chain array has much higher degrees of freedom of 
the ground states and low-lying excitations than those of 
the honeycomb model. This makes the difficulty to peel 
out a given set of the anyons [l^. Carefully checking 
the mapping, one sees that this is caused by the lost of 
the couplings between the chains when a Jordan- Wigner 
transformation is used in the mapping. The source of 
the lost of the coupling arises possibly from the unproper 
treatment of the periodic boundary condition in Jordan- 
Wigner transformation. An additional boundary Hamil- 
tonian should be added [l^, which is a bit complicated. 
Focussing only on the low-lying sector, we can introduce 
vertical Ising ferromagnetic couplings in two given ver- 
tical chains (the 'boundary chains') to approximate the 
boundary couplings. After modifying the totally decou- 
pled chain array model to the boundary coupled chain 
array model, the degeneracy of the ground states be- 
comes consistent with the honeycomb model and so are 
the low-lyings excitations. Therefore, this boundary cou- 
pled Ising chain array model is equivalent to the toric 
honeycomb model in the low-lying sector. 
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FIG. 1: (Color online) Toric honeycomb lattice. The green 
loop (armchair) and the red loop (zig-zag) are two indepen- 
dent non-trivial loops. 
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II. TORIC HONEYCOMB MODEL 

The model we propose is very simple. Consider a hon- 
eycomb lattice in a torus (See Fig. [T]). The number of 
total sites is 2 A'' and each site is occupied by a 1/2-spin. 
Define two plaquette operators 

Wp = afa^alala^.a^ 

Wp = (T^a|a|a|a|a6^ (f) 

where 1, • • • , 5, and 6 are the site indices of a plaque- 
tte P as shown in Fig. [T] It is easy to directly check 
Wp = Wp = 1, i.e, the eigenvalues of them are ±1; and 
[Wp,Wp/] = for all plaquette. The Hamiltonian is 
simply defined by sum over all plaquette 

p p 

where Ji > and J2 > and thus Wp and Wp' are 
two-independent integrals of motion. This spin model 
is invariant in xy z permutation. This model is exactly 
soluble: the ground states are given by the states with 
Wp = Wp = 1 for all plaquette and the excited states 
are of at one plaquette having Wp = — 1 or Wp = — 1, 
which will be studied later in details. 



III. GROUND STATES AND DEGENERACY 

The ground states are of the form 

\G) ^Y[il + Wp)il + Wp)\<f>), (3) 
p 

where \(j>) is an arbitrary reference state, e.g., |(/)) = | t 
• • • t) with each 'f' denoting the eigenvalue 1 of cf . The 
number of \(f>) is 2^^, the dimensions of the whole Hilbert 
space . The above construction of the ground states 
means the ground state space L = J^/G where G is 
the group generated by all independent Wp and Wp. 
On the torus, due to the periodic boundary condition, 
there arc two constraints Y[p Wp = Y\p Wp = 1 for 
the generators, which means G ^ (Z2)2^-2. Therefore, 
according to the general theory given by Gottesman (2^ , 
the degeneracy of the ground states are 2'^^~^^^~^^ = 4, 
which is exactly the same as that of the Kitaev toric code 
model 

One can even make a more intuitional proof similar 
to what Kitaev did For a plaquette P, one de- 

fines Kij to be a-a^ for = (1,2), (4, 5), crfcrj for 

(i,j) = (2,3), (5,6), and afa^ for (i,j) = (3,4), (6,1); 

to be afa| for (z, j) = (1,2), (4,5), afa] for = 
(2, 3), (5, 6), and afaj for (*,j) = (3, 4), (6,1). For any 
closed loop, one can define two different operators 

C = K^jKjk ■ ■ ■ K,^, C = K^JKJk ■ ■ ■ Ksf (4) 



For any trivial loop on the torus, C and C are iden- 
tical to the identity. There are only two unequivalent 
non-trivial loops as shown by the green and red in Fig. 
[H Along the green loop, Cg = a^a^a^afa'^a^a^a^ ■ ■ ■ 

and Cg = af aj a^af af^a^a^af ■ ■ ■ while along the red 
loop, they are Cr = crfcrjcr^afa^a^a^af ■ ■ ■ and Cr = 

'^f'^j'^k'^f'^m'^n'^s'^t " ' ■ These four operators generates 
all linear operators acting on the ground state space. In 
this way, we explicitly prove that the ground states are 
fourfold degenerate. 

IV. LOW-LYING EXCITATIONS 

A. Excitations 

Low-lying excitations of this exactly soluble model are 
given by 

<\G), 

n 41^)' 

i'g<ib 

Xn\G)=al n '^t'JG), 

<7^P^\G)^alal_,al_,---\G), 
4S'|G)=«-i<-3---|G), 
4S'|G)=«-i<-3---|G), 
4V|G)=af^_2<_4<_6---|G), (5) 

where the sites it, are shown in Fig. [T]and P [P') label 
the plaquette where the site T'=i;, (if, — 2); and the order 
of the sites is defined as follows: is > jt if the zig-zag line 
including is is higher than that of jt or if is is on the right 
hand of jt when they are in the same line. Except the 
first local excitation, all others are string-like excitations 
which are plotted in FiglH 

Notice that ct^, cr^ and anticommute with each 
other and therefore, if there is a red triangle (cr^) or 
black triangle (cr^) at site 3 or site 6 in a plaquette P, the 
eigenvalues of Wp and Wp of this plaquette become — 1 ; 
if there is a blue circle (cr^) at site 1, site 2, site 4 or site 
5 in P, the eigenvalues of Wp and Wp of this plaquette 
also become —1; if there is a red triangle at site 1 or site 
4, or a black triangle at site 2 or site 5, the eigenvalue 
of Wp becomes —1; if there is a black triangle at site 1 
or site 4, or a red triangle at site 2 or site 5, the eigen- 
value of Wp becomes —1. Taking ipi^ as an example, the 
blue circle at site 1 in plaquette A (See Figl2]) makes the 
eigenvalues of Wa and Wa become — 1 , while the red tri- 
angle at site 5 makes the eigenvalue of Wa become — 1, 
and finally, the eigenvalue of Wa in plaquette A becomes 
— 1. For plaquette B, the red triangle at site 1 makes the 
eigenvalue of Wb become —1. It is easy to check that 
the eigenvalues of W and W in other plaquette are not 
changed. 
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FIG. 2: (Color online) The low-lying excitations. From left 
to right and up to down, they are ipi^ , , (jp ' , Cp p, , o"p p/ , 

and (jpV. The up-triangles (red), down-triangles (black) and 
blue circles indicate acting the operators , o'" and at 
these sites on the ground state, respectively. 

Although string-like excitations look like non-local , 
their energies are located at one or two sites. The energies 
of excitations can be calculated directly, which in turn are 
4Ji, 4J2, 4J1+4J2 for V'.X,^^ and 2J1+2J2 for ct(1'2,3)^ 
respectively. 

Notice that changing the subscript i}, to i^, the exci- 
tations do not change, i.e., we can redefine excitations as 
follows 

aW|G)=«_2<-4---|G), 
4%,|G)=«_i<,_3---|G), 
4%'|G>=«_i<,_3---|G), 
a(^)|G)=<_2<-4<-6---|G>. (6) 

The only difference is the interchange between and 
(7^, and it is easy to check excitations defined in Eqs. ([5]) 
are the same as those in Eqs. 



B. Fusion Rules 

All excitations are explicitly expressed by Pauli ma- 
trices. The fusion rules of these operators then can be 
directly calculated (See TABLE|l|. 



TABLE I: Fusion rules of excitations. 
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The fusion rules of the closed subset 
{I,ipi^,ap\ap^p,} are exactly the same as the fu- 
sion rules in Kitaev toric code model if we identify 
7/), CT^^^ and cr^^^ as e, e and m[l|, [13] : 

^2^(a(l))2 = (a(2))2 = l, 

(7(1^(2)=^-, =a(2), a(2V = a(i). (7) 

Here '= ' may be up to a sign and/or 'i\ Similarly, 
the subset {/, Xh , CTp ^ , Cp^p, } also obeys the same fusion 
rules. 

Because of the periodic boundary condition for a torus, 
the number of string-like excitations must be even. 

V. BRAIDING RULES AND ANYONS 

Enlightened by this equivalence in the fusion rules, we 
check the braiding matrix. At first we find the statistics 
of the same kind excitations. It is easy to get [dp -* , ap)] = 
0) [^Pi'',p' ' CTpj.pJ = 0, and {4>ib,4'jb} = 2(5^, i.e. Cp' and 
Op'pi obey Bose statistics, while V'it obey Fermi statis- 
tics. Since tr^^V^^^ = ')/', the interchange between two '0, 
i?^^, can be considered as pro duct of R^{\)^{\) , R„(2)„{2) , 
R^{i)^{2) and R^(2)^(i), i.e.,[2lj 

R^ll> = -Ro-(l)(T(l)-R(T(2)(T(2)-Rcr(l)cr(2)-RCT(2)o-(l)- (8) 

Clearly, R.^^ = -1, i?^(i)^(i) = 1 and i?^(2)^(2) = 1, so 
Rg.{i) ^(2) R^(2) ^{1) = —1. This means when ct^^-* circles 
around cr'-^^ or vice versa, a minus sign is acquired. In 
other words, they obey mutual semionic statistics. Braid- 
ing ipii^ with cr^^^ or cr^^-' also gives —1. Similarly, statis- 
tics among XibJ c''*'^'' and f^^'' is also semionic. 

VI. MAPPING TO A SQUARE LATTICE 

A. Jordan- Wigner Transformation 

This model can be mapped to a square lattice via a 
Jordan- Wigner transformation. For simplicity, we first 
restrict to the open boundary condition. We have defined 
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4'ib-i Xibi and Xiu, before. It is easy to check they 
are Hermitian and {ipi^^^tpj^} = 25ij, {V'i„)V'j„} = 25^ 
and {ipitT^j^} — 0. These relations are valid for x '"^ 
well. Furthermore, ■0 and x are anticommutative. These 
operators can be thought as Majorana fermions, and can 
be combined into 'complex' fermions 



1 



1 



1 



1 



(9) 



Expressing Pauli matrices by ip and X: we get 

= n '^^'s ' = ^^i" n ' 
<=x.„ n <=V'H n <■ (10) 

Thus, Wp and VFp can be expressed by i}} and x as 

Wp ^iXt.XzJXnXj^- (11) 
Using the complex fermions, one has 

Wp = (2nt,,; - l)(2nt,, - 1), 

M^P = (2n;,, - l)(2n;,, - 1), (12) 

where Ug^i = cl .^Cg^i are the fermion number operators. 
The Hamiltion can be rewritten as 

H^-Ji (2"T,^ - l)(2%j - 1) 

{ij)hd 

-h J2 (2"^.' - l)(2"ij - 1)- (13) 

{ij)hd 

The symbol {ij)hd means the sum is over nearest neigh- 
bors along the horizontal diagonals of squares (See FigEl). 
In this way, we transfer the honeycomb lattice model to 
square lattice model composed of series of decoupled Ising 
chains as shown in FiglSl 



B. Ground states and Excitations 

This square lattice model has been studied in our pre- 
vious work [l^ . The ground states of this Hamiltonian 
are 2^^-fold degenerate, i.e., every individual chain is fer- 
romagnetic, i.e., for a set of spins {si, • • ■ , Sa, • • • , Sm}, 



M 



iGM)=niG^»)=n^i.jo)' 



(14) 



a=l 



where M is the number of chains, ia is the site index in 
the a-th chain. The low-lying excitations above a given 




Ib^Jb 



FIG. 3: (Color online) Map from the honeycomb lattice to 
the square lattice. Ferromagnetic couplings are added to the 
two red lines. 



ground state |G{s}), for a given site ia, reads 



IG 



^p\G{s}) = n ^^'JG{s})= n *^«..^'.4..JGt,}), 

b — " b — " 

W^,p,|G|,}) = n J-,.H,JG{,}) 

= n *c^->^'l4„«^c3„,,JG{,}), (15) 



ii<ia 



where P and P' denote two plaquette on the right and 
left of ia, respectively, s =1 (t) if s =t (i)- The order 
of sites is defined by < ia if i'l, is on the left hand 
of id for 6 = a or i'j, is in a chain lower than the chain 
with ia- 'H,'D,J- create a hole, a double occupant, and 
a spin-flip. W, and W'* create a half-infinite string 
of spin-flips, a spin-flip string with a double occupant 
and a spin-flip string with a hole, respectively, since the 
spins of fermions at sites ij, < ia are flipped from their 
ground state configuration while those at jc > ia keep in 
their ground state configuration. These excitations are 
equivalent to those on honeycomb lattices with 



C. Equivalence in Low- lying Sector 

As we have seen, the degeneracy of the ground states 
of the decoupled Ising chains is much higher than the 
toric honeycomb model. This means two models are 
not equivalent. This is caused by unproper treatment 
of the periodic boundary condition of the toric lattice 
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in the Jordan- Wigncr transformation. Due to the peri- 
odic boundary condition, an additional boundary Hamil- 
tonian Hb should be added to Eq. (fTS)) [l9| . The exact 
form of Hb is complicated. Instead, we can use an ap- 
proximation to deal with the periodic boundary condition 
so that the low- lying excitation sector is correct. We in- 
troduce Ising ferromagnetic boundary couplings to the 
system, i.e., for the two vertical red lines in FigjSl we ap- 
proximate the boundary effect from the Jordan- Wigner 
transformation through a boundary coupling Hamilto- 
nian 

Hb = -gi ^ n.isnjs - 32 ^ nisHjs , (16) 

where 1 and 2 label two vertical red lines in Figl3] and 
{ij)i,2 mean the sum is over corresponding vertical lines. 
We assume gi^2 is much larger than J 1^2, the energy scale 
of the low- lying excitations, so that the opposite nearest 
neighbor spins along a given boundary line do not cre- 
ate additional low-lying excitations. Clearly, this kind of 
couplings make all odd/even chains have the same states 
in a given ground state. Therefore, the degeneracy of the 
ground states is the same as two decoupled Ising chains, 
i.e., 4. And couplings will not change the ground state. 



Due to the periodic boundary condition, low-lying exci- 
tations will appear in pairs. Thus, as long as excitations 
do not cross the vertical boundary lines, excitations are 
low-lying . The statistics of excitations in square lattice 
are the same as those in honeycomb lattice, see p^ . 



VII. CONCLUSIONS 

We have constructed a toric honeycomb model and 
proved that it is equivalent to Kitaev toric code model 
in square lattice. This honeycomb model is helpful to 
understand the toric code of anyons in a more symmetric 
way. It was also shown that this model can be mapped 
into an Ising chain array with boundary couplings, which 
proposed a possible realization of the anyon models in a 
conventional interacting fermion model. 
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